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, Q [ Abstract. We develop an analytical tool which is adept for detecting 

. shapes of oscillatory functions, is useful in decomposing homogenization 

problems into limit-problems for kinetic equations, and provides an effi- 
, cient framework for the validation of multi-scale asymptotic expansions. 

We apply it first to a hyperbolic homogenization problem and transform 
it to a hyperbolic limit problem for a kinetic equation. We establish con- 
ditions determining an effective equation and counterexamples for the 



PL, . , 

■ case that such conditions fail. Second, when the kinetic decomposition is 



applied to the problem of enhanced diffusion, it leads to a diffusive limit 
problem for a kinetic equation that in turn yields the effective equation 
of enhanced diffusion. 
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>■ ' 1. Introduction 

in ' 

CN ■ Homogenization problems appear in various contexts of science and engi- 

neering and involve the interaction of two or more oscillatory scales. In this 
work we focus on the simplest possible mathematical paradigms of periodic 
homogenization. Our objective is to develop an analytical tool that is capa- 

■ ble of understanding the shapes of periodic oscillatory functions when the 
scales of oscillations are a-priori known (or expected), and use it in order 
to transform the homogenization problem into a limit problem for a kinetic 
equation. The calculation of an effective equation becomes then an issue 

^ ■ of studying a hyperbolic (or diffusive) limit for the kinetic equation. The 

■ procedure is well adapted in identifying the specific characteristics of the 
underlying homogenization problem and provides an efficient tool for the 
rigorous justification of multiscale asymptotic expansions. 

The main idea is motivated from considerations of kinetic theory. When 
the statistics of interacting particles is studied it is customary to introduce 
an empirical measure and to study its statistical properties in the (weak) 
limit when the number of particles gets large. Likewise, for an oscillating 
family of functions {u £ } if we want to study the shape of periodic oscillations 
at a predetermined scale we may introduce an inner variable that counts the 
content of oscillation at such scale. For instance, to count oscillations at the 
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scale I one can introduce 

(1.1) f £ (x, v) = U £ (x)5 p (v - ~) 

where 5 P is the periodic delta function, and study the family {/ e }. A-priori 
bounds for {u £ } translate to uniform bounds for {/ e }: if for example u £ is 
uniformly bounded in L 2 , u £ L 2 , then f £ Gj, L 2 (M(T d )) and, along a 
subsequence, 

(1.2) f £ ->> / weak* in L 2 (M(T d )) , 

where M(T d ) stands for the periodic measures. In addition, the resulting / 
is better: / G L 2 (L 2 (T d )). 

The above object should be compared to the concept of double-scale limit 
introduced in the influential work of Nguetseng [T7j and applied to a variety 
of homogenization problems [U [151 M] ■ In the double-scale limit one tests 
the family {u £ } against oscillating test functions and develops a represen- 
tation theory for the resulting weak-limits. It turns out, [17], that for a 
uniformly bounded family u £ G& L 2 and test functions ip periodic in v 

(1.3) J u £ (x)(p(x, — )dx — > y y f(x,v)(p(x,v)dxdv 

where / G L 2 (L 2 (T d )). The reader should note that this is precisely the 
content of (jl.ll) . (|1.2I) . which thus provide an alternative interpretation to 
the double scale limit. However, what seems to have been missed, perhaps 
because Nguetseng's analysis [IT] proceeds without writing down (11. lj) but 
rather by establishing directly (jl.3p . is that the measures f £ satisfy in their 
own right very interesting equations. This is a consequence of additional 
properties, like 

(1.4) (v x + -V v )fe{x,v) = V x u £ (x)5 p (v - - £ ) , 

obtained by applying differential operators that annihilate the singular mea- 
sure. Properties like (jl.4p . in turn, suggest a procedure for embedding ho- 
mogenization problems into limit problems for kinetic equations. In the 
sequel we develop this perspective, using as paradigms the problem of hy- 
perbolic homogenization, and the problem of enhanced diffusion. 

The double-scale limit |17j along with the technique of multiscale asymp- 
totic expansions [5] have been quite effective in the development of homoge- 
nization theory with considerable progress in several contexts {e.g. [16], [T], 
0]) [2], [IH]) [IS]). Other tools have also been used for the homogenization 
of linear hyperbolic problems: Among them are of course Young measures, 
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developed by Tartar and used for the homogenization of some particular 

linear transport equations in two dimensions (see [20] and [21] ). Wigner 

measures (see [IT]) may also be mentioned. 

As our first example we consider the hyperbolic homogenization problem 

du £ / Xs. 

— -+a{x, -) -VxUe = 

(1.5) dt e 

u £ (0,x) = U°(x,~), 

with a(x, v) a divergence free field periodic in v, is transformed to the prob- 
lem of identifying the hyperbolic limit e — ► of the kinetic initial-value 
problem 

Qf 1 

— -f + a(x, v) ■ V x f £ + -a(x, v) • V„/ e = 0, 

(1-6) dt 

f £ (t = 0,x,v) = U°(x,v) 5 p (v--) 

Homogenization for (|1.5|) has been studied by Brenier [6] , E [9] , Hou and Xin 
[15] and, in fact, the effective equation is sought - motivated by the double- 
scale limit - in a class of kinetic equations. Eq. (|1.5|) is by no means the only 
interesting hyperbolic problem for homogenization; we refer to [3], [13], [12] 
(where a kinetic equation itself is homogenized), and to [2] for an example 
concerning a Schrodinger equation (the list is of course not exhaustive) . 

For (II. 5j) . our analysis proceeds by studying the hyperbolic limit for the 
kinetic equation (jl.6p . We find that if the kernel of the cell-problem 

(1.7) K x = |g G L 2 (R d x T d ) a(x,v)-V v g = 0inV'\ 

is independent of x, then it is possible to identify the effective equation. 
Namely, when the vector fielfd a = a(v) is independent of x the effective 
equation for / reads 

?r + ( Pa ) • = 

(1.8) dt y ' J 

f(t = 0,x,v)=PU°(x,v), 

where P is the projection operator on the kernel K, and in turn u = L, d fdv 
(see Theorem 13. ip . By contrast, when a = a(x) and K x depends on x, a 
counterexample is constructed that shows that the effective equation can not 
be a pure transport equation (see section [3T2|) . In section this analysis is 
extended for homogenization problems where a periodic fine-scale structure 
is transported by a divergence-free vector field (see equations (|4.ip and (|4.4p ) 
analogous results to the case of (11.51) are found. Such kinetic equations might 
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turn very useful for devising computational algorithms for the computation 
of homogenization problems. 

A second paradigm is the problem of enhanced diffusion 

1 x 

d t u £ + -a(x, -) • V x u £ = aA x u e 
(1.9) e e 

u £ (0,x) = U°(x,^) 

with a(x, v) periodic, divergence-free and with mean Jj d a = 0. The results 
formally obtained by multiscale asymptotics have been validated for this 
problem by McLaughlin, Papanicolaou and Pironneau [16) . Avellaneda and 
Majda [4J, and Fannjiang and Papanicolaou [10) . We revisit this problem 
from the perspective of the kinetic decomposition and transform it to the 
problem of identifying the e — > limit 

+ -a(x, v) • V x f £ + v) • V v f £ - Q>A v f e ) 

( L1 °) = aA x f £ + — V x -V v f £ , 

e 

f e (t = 0,x,v) =U°(x,v) S p (v-^). 

The latter is a limit for the transport-diffusion equation (jl.lOp in the so- 
called diffusive scale, and its analysis provides the effective equation (|5.7[) - 
(|5.9|) of enhanced diffusion (see Theorem 15.1)) . This example indicates the 
efficiency of this approach in the rigorous validation of multi-scale asymp- 
totic expansions. 

Finally, we note that the scales of the drift and of the diffusion in (|1.9p 
may be chosen differently from 1/e and a, yielding other interesting homog- 
enization problems, see for instance Capdeboscq (7) [8]. 

The article is organized as follows. Analytical considerations like the 
proper definition of (jl.ll) , the characterization of the weak limit points of f £ 
under various uniform bounds, the differential relations such as (jl.4p . and 
the identification of asymptotics for f £ are developed in section [2] and in 
appendix I. In section [3) we study the hyperbolic homogenization problem 
(|1.5p . derive the effective equation, and produce the counterexample men- 
tioned before. Some material from ergodic theory needed in the derivation 
is outlined in the appendix II. In section HI we study the transport via a 
divergence-free field depending on an oscillating fine-scale, we derive the 
associated kinetic equation, and discuss the connection of the two formula- 
tions via characteristics and the derivation of an effective equation. Finally, 
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in section [5] we study the parabolic homogenization problem (|1.9p and derive 
the enhanced diffusion equation via the kinetic decomposition. 

2. Multi-scale decomposition 

Let {u £ (x)} be a family of functions defined on a open set Cl C M. d that 
contains periodic oscillations and suppose that the scales of oscillations are 
either a-priori known (or anticipated). Our goal is to introduce an analytical 
object that will prompt the anticipated scale(s) of oscillations and quantify 
the structure of oscillations in the family at the preselected scale(s). 

Suppose that periodic oscillations of length e are anticipated in the family 
{u £ }. To focus on them we consider a periodic grid with sides of length e in 
each coordinate direction. The grid splits the Euclidean space into distinct 
cubic cells of volume £ d , and it is arranged so that the centers of the cells 
occupy the lattice eL d . Let Vt be placed on that grid, and define a function 
Xe '■ CI —> ^ d that maps the generic i £ !1 to the center Xe{ x ) °f t ne 
cell containing x. To each point x € Cl there is associated a decomposition 
(Xe{x),v) where Xs( x ) £ eL d stands for the center of the cell that x occupies, 
and v G T d is the vector difference x — Xe{ x ) as measured in units of distance 
e, that is x = Xe( x ) + sv. We introduce the quantity 

(2.1) f £ {x,v)=u £ {x)5 p {v- X ~ X ^ x) ), x£R d ,v£T d , 

where 5 P stands for a periodization of the usual delta function with period 
1 in each coordinate direction, and T d stands for the d-dimensional torus, 
the quotient of M. d by the subgroup Z . 

We note that the map x \— * (xe( x ), v ) is single valued for points that fall 
into a single cell, but multi-valued for points that fall onto the boundaries 
between adjacent cells. For the latter points there would be two different 
decompositions (xe, v ) an d {x' s -,v') associated to the same point x £ Q. 
Nevertheless, in that case x = \e + £v = x' £ + £v ' an dj due to the use of a 
periodic delta function, 

Hence, both decompositions provide the same outcome in (|2.ip with f £ de- 
fined for and v £ T d . 



6 



P.-E. JABIN AND A.E. TZAVARAS 



The operator V x + -\/ v annihilates the form v — | and that - at least 
formally - yields the formula 

(2.2) (y x + l Vv )f £ = (y xU£ ) Sp (v-~) 

In the sequel, we provide formal definitions for the decomposition (|2.ip and 
extensions as well as differentiation properties like (|2.2p that are helpful in 
later sections for validating multiscale expansions. 

2.1. Definitions. We make extensive use of distributions defined on the 
torus T d . Such distributions are in one-to-one correspondence with periodic 
distributions T on of period 1 in each coordinate direction, that is dis- 
tributions satisfying for i = G Z d the property t{F = T where 
Tj is the shift operator, see [El p. 229]. The same notation is used for 
both interpretations of periodic distributions. Let 5 P be the periodic delta 
function of period 1, defined by its action < 5 p ,ip >= ip(0) on continuous 
periodic test functions ifi £ C(T d ). 

We use the notation C p = C(T d ) for the continuous periodic functions, 
C oc _ c°°(J d ) for periodic test functions and M p = M 1 (T a! ) for the periodic 
measures, with period 1 in each coordinate direction. Recall that C p is 
separable and that bounded sets in M p = (Cp)* are sequentially precompact 
in the weak-* topology of M p . 

2.1.1. The double-scale kinetic decomposition. Our objective is to define the 
product 

(2.3) f e (x,v) = u £ (x)5 p (v - -) 

which, in view of the periodicity of 5 P and ^ G Z d , coincides with (|2.ip . 
Products of distributions are not in general well defined. However, due 
to the tensor-product-like structure the product in (|2.3p is well defined by 
using the Schwartz kernel theorem |14[ Ch V]. We defer the details for the 
example in section 12.1.21 an d note that the definition of (j2.3|) is effected for 
u £ G L} oc (p) by setting 

(2.4) <f £ ,0>= [ u e (x)6(x,-)dx 

for 6(x,v) G C£°(£l; C°°(T d )). Moreover, we have the differentiation for- 
mula: 
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Lemma 2.1. For u £ G 

j / X \ X 

(V x + -V„) [u £ 5 p (v - -)) = (V x u £ )5 p (v - -) 
Proof. For G C c °°(0; C£°), we have 

< (Va + -v v )f E ,e> = -<f £ , (v x + ^v„)# > 

= - / u E (x) (v x e + -v v e) (x, -) dx 

jRd, \ e / e 

u £ (x)V x .(6>(x,^))<ix 

(V x n £ )(x)6*(x, -) 
< (V x u £ )5p(z;- -),e > 



□ 



2.1.2. A generalization. Let a(a;) be a smooth vector field and u E G L\ oc {Q). 
We proceed to define the product 

(2.5) f £ = U£{x)6p ( v -^l). 

Naturally it should act on tensor products (p (g> ip of test functions via the 
formula 

(2.6) <f £ ,y®xjj>= f u £ (x)<p(x)iP(^-)dx. 

JM. d e 

To define (|2.5p . we employ the Schwartz kernel theorem [14| Thm 5.2.1]. 
Consider the linear map 

K : C°°{T d ) -^V'(Q) defined by Kip = u £ (x)ip(^-) 

If tp n — > in C°°(T d ) then /C?/Vi in 2?'(J7). The kernel theorem implies 
that there exists a unique distribution K such that < K,<p®^p >= (ICip, p), 
that is K acts on tensor products via (|2.6p and is the desired product. It 
satisfies, for 6 G C c °°(ft; C7°°(T d )), 

(2.7) <UeSp (v-^),e>= [ u £ (x)e(x,^)dx, 

£ JR d e 

which can also serve as a direct definition of / £ . Of course smoothness of a(x) 
is required for the above definition: at least a G C(J7; U. d ) if f £ is interpreted 
as a measure, and more if f e is interpreted as a distribution and we need to 
take derivatives. 
We now prove. 
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Lemma 2.2. Let u £ £ Wj£(Q.) and a £ C 1 ^; T i ). Then 

(V« + - £ (Va) T V v ) Ud p (v - ^)) = (V x u e )6 p (v - ^1] 
Proof. For the k-th. coordinate, we have 

<(^£^0«*<.-^).»> 

= -<«A(»--7 i )'V + iE». J (^)> 

j 

6*(x, ) ) 

£ / 

(d x . u e )(x) 0(x,^-)dx 
i e 



□ 



2.1.3. A multiscale kinetic decomposition. We pursue next the construction 
of decompositions in cases when more than two scales are involved. Suppose 
that for an oscillating family {u £ } we wish to focus on oscillations at the 
scales 1, | and p-. We define 

(2.8) f £ (x,v,w)=u E (x)6 p (v-~)6 p (w-^), x £ fi, v £ T d , w G T d , 

or, in terms of the action on test functions, 6(x, v, to) £ C^°(Q; C°°(T d x T d )) 
via the formula 



(2.9) < 



f £ ,0>= / e( X ,-,^)dx. 



In a straightforward generalization of Lemma 12.14 fe satisfies, for u e £ 



Wj j}(Sl), the differentiation formula 



(2.10) (V,. + ^V,, + ^V w )f e = (V x u £ ){x)5 p (v - ^)5 p (w - ~) . 

To motivate the definition (|2.8p consider for simplicity the case that 1/e 
is an integer. Fix a first grid of size e and introduce the quantities Xei x ) 
and v = x ~ Xe £ T d as before. To focus on the scale \ we consider a second 
grid of length e 2 embedded in the first grid. When 1/e is an integer, the 
grids fit perfectly onto one-another. Define the function ip £ : T d — ► eZ d 
that takes the generic point t> to the center of the inner cell containing v, 
and introduce a second inner variable to = - ^ e — £ T d describing the 
vector distance between v and the center of the inner cell containing v in 
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Ue(x)Sp(v- -)6 p (w- -) , x G ft, v G T d ,w G T d , 



units of length e. The process defines a decomposition of the physical space 
x i — > (xe(x), v, ip e (v),w), and allows to define a kinetic function representing 
three scales by 

f £ (x,v,w) = u £ (x)5 p (v - £ X£ )S p (w - - 

This definition is also good when 1/e is not an integer as can be seen by the 
formula (|6,2|) in the appendix. 

2.2. Multiscale analysis of uniformly bounded families of functions. 

Nguetseng |17] introduced the notion of double scale limit, which has been a 
very effective technical tool in the development of periodic homogenization 
theory. His approach does not use the kinetic decomposition (12. ip . but the 
double-scale limit is precisely the weak limit of the measures introduced in 
(|2.ip . For this reason, we will review the results of Nguetseng [17] from 
the perspective of the theory presented here, and produce some further 
asymptotic analysis of kinetic decompositions for uniformly bounded families 
of functions. In the sequel, the notation u £ G& X means that the family {u £ } 
belongs in a bounded set of the Banach space X. 

2.2.1. Uniform I? '-bounds. Suppose first that {u £ } satisfies u £ G^ L 2 (S7). 
We define f £ by (|2.3p and note that 

(2.11) f £ e b L\n;M p ). 

The Riesz representation theory asserts that there is an isometric isomorhism 
between the dual of C p = C(T d ) and the Banach space of periodic Radon 
measures M p = M l (Y d ) on the torus. Since C p is separable, bounded sets in 
Mp are sequentially precompact in the weak-* topology of M p . Also, since 
C p is separable, so is L 2 (fL;C p ) and thus bounded sets in L 2 (Q;M p ) are 
sequentially precompact in the weak-* topology of L 2 (f2; M } 



pi- 



As a consequence (]2.1ip implies that, along a subsequence, 

(2.12) f E ->> / weak-* in L 2 (VL; M p ) 

with / G L 2 (Q; M p ), that is 
(2.13) 

< f £ ,0> = J u £ (x)6(x, 2) dx 

/, 6 >= If f(x, v)6(x, v) dxdv for 9 G L 2 (ft; C p ) 
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Examples. A few examples will illustrate the properties of this convergence. 

1. Note first that 

(2.14) 5 p (v-^) ^1 weak-* in L°°(ft; M p ) , 
that is 

(2.15) < S p (v- -),6 >= [ 9(x,-)dx -> f [ 6(x,v)dxdv, 

for 9 £ L 1 (0; C p ). This is a classical result, see [5], and a proof is provided 
for completeness in the appendix. 

2. If u £ — ► u strongly in L 2 (f2), then 

(2.16) u e (x)5 p (v - ~) u(x) weak-* in L 2 (ft; M p ). 

Indeed, since ii e — > u in L 2 (f2) and ^(f ) — J T d ip(v)dv weakly in L 2 (f2), for 
9 = (p ® if) a tensor product 

oc fx 
< u e 6 p (y ),<p<S)ip> = / u £ (x)^(— )(f(x)dx 



u{x)ip(x)ij)(v)dxdv 

Finite sums of tensor products ^ ■ (pj®tpj are dense in L 2 (Q; C p ) and (|2.16p 
follows by a density argument. 

3. For u £ = a(~ ). where a(t>) is a periodic function, we obtain 

(2.17) oC-)-" f a(v)dv 
and, by (|6.ip . 

(2.18) a (£)«j p ( w _£) _. a ( w ) 

Observe that the weak limit (12. 17ft retains only the information of the aver- 
age of a while the double scale kinetic limit (|2.18p also retains the informa- 
tion of the shape of a. Equation (|6.2p in the appendix indicates that if we 
were to prompt u £ = a(|) with test functions oscillating on a scale different 
than e then the information on the shape of a is lost and only the average is 
perceived. Therefore, the double scale decomposition works well when the 
scale of oscillations are a-priori known, so that the right oscillatory scale is 
prompted. 
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4. Finally, let u £ = a(~ )£>(;§■) where a and b are periodic functions. Then 
EE]) and (jfT2|) give 



X X f f 

a (-) & (^) / a (y)dy / 

■Z/ Xj Xj i 

a {-)K—)bp{ v ) a ( v ) / b(z)dz 



e e' J y<i 

a(~)b(^)5 p (v - ^)6 p (w - -) a(v) b(w) 



In the following proposition, we give a simplified proof of 17, Thm 1] 
concerning the double-scale limit for sequences that are uniformly bounded 
in L 2 {n). 

Proposition 2.3. Let u e G& L 2 (S7). Then, along a subsequence, 
fe—^f weak-* in L 2 (f2 ; M p ) 

with f G L 2 (!Jx T d ). 

Proof. Let 9 G C c °°(ft; C p °°) be a test function. Then 

< f £ ,0>= / u £ (x)6(x, -)dx 
Jn e 



and 



< / £ ,0 > | < ||u E 



Ii2(n) (J^\0(x,^)\ 2 dx 
<C(J \d(x,^)\ 2 dx) 



Hence, / £ G b (L 2 (ft; C p ))*, / £ / weak-* in L 2 (Q; M p ) and / G L 2 (ft; M p ). 
Moreover, 

,Ll /)g>l =lim!</ e ' g> '<C 



lFllL 2 (nxT d ) lrllL 2 (nxT d ) 
and / G L 2 (ft x T d ). □ 

2.2.2. Uniform H 1 -bounds. Next consider the case of families {u £ } that are 
uniformly bounded in H 1 ^). The first proposition is essentially a rephras- 
ing of [13 Thm 3]. 
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Proposition 2.4. Let u £ G& H l {Q). Then, there exist u G H (0), 7r G 
L 2 (J7; ^T 1 (T d )) smc/i i/iot, along a subsequence, 

f £ = u e 8p(y ) — v u(x) weak-* in L 2 (S7; M p ) 

(V x u e — X7 x u)5 p (v ) — 1 V„7r(:z;,f) weak-* in L 2 (£l;M p ) 



Proof. Along subsequences (whenever necessary) u e — it weakly in 

u e — > n in L 2 (f2) and u £ 5 p (v — — 1 u weak-* in L 2 (f2;M p ). Moreover, 

Proposition 12.31 implies 



JW - -J 9i{z,v) 



dxi dxi 

weak-* in L 2 (n-,M p ) with # G L 2 (0 x T d ). 
For y> G C c °°(0), ^eC™, we have 



(u e - u) 



dxdv 



(p(x)ipi(-) dxdv 
dxi v s y 

<7i(x, v)ip(x)vjji(v)dxdv 



We apply the above formula to a test function ^ = • Vv) that satisfies 
div„* = 0. Then 



(u £ - u) 



(fxcfo — > 



and we conclude that for a.e. 

(2.19) / \g i {x > v)ij} i {v)dv = Q for any ^ with div„\I/ = 0. 

A lemma from [17\ Lemma 4] then implies there exists tt G L 2 (Q; H 1 (T d )) 
such that G = (g±, ...,<?d) = V„7r. □ 

The next proposition is novel and establishes the asymptotics of f £ when 
the family {u £ } is uniformly bounded in H 1 ^). 

Proposition 2.5. Let u £ G& H l (£i). Then 



(2.20) 



^(6 p (v - ~) - 1) ^ inV 
g £ := u £ {x)- e {5 p (v - ^) - 1) G 6 H^^Mp) 



g £ — 1 5 weak-* in H M p ) 
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and f £ enjoys the asymptotic expansion 

f e = u £ + eg + o(e) in V . 

Proof. It is instructive to first give a quick proof for the case of one space 
dimension. Consider the function H{v) = v , v S [0, 1] and let H p denote its 
periodic extension of period 1. H p £ L°°(R) satisfies d v H p (v) = 1 — 8 p {v) 
and 

(2.21) l( Sp (v-~)-l)=d x H p (v-~) 
Using standard properties of weak convergence we obtain for 9 G C^(Q, C p ) 
< H p (v ),#>—► / / / H p (v — w)9(x,v)dwdvdx 

£ ./<> J1 d JT d 

and 

< - £ {t P {v-- e )-i),e>=<d x H p {v-- e ),e> 

= - <H p ( v -^),d x e> 

— v — / / / H p (v — w)d x 9(x,v)dwdvdx 
Jn Jj d Jj d 

= 0. 

Using lemma I2TT1 and (|2.21|) we have, for u e E the identities 

1 X X 

(d x + ~d v ) (u £ H p (v - -)) = (d x u £ )H p (v - -) 

and 

u £ -(5p(w ) - l) = d x (u £ H p (v ) - (d x u £ )H p (v ) . 

£ £ V £ / £ 

The remaining three properties are direct consequences of the last formula. 

Consider next the multi-dimensional case. For 9 £ C^(fl;C p ), let supp# 
denote the support (in x) of 9, and fix e < -i dist (supp 0, 50). We cover 
supp# by cubes centered at points Xk G ^ d of latteral size e. The 
number of the cubes covering supp# is of the order e _d O(|supp 9\), and the 
covering is arranged so that supp# C U^ =1 Ck C 0. We observe that by 
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construction jXk £ and compute 



< g £ ,6 >= / u e (x)- (#(iE, — ) — / f cfe 
Jn e \ e J T d ) 

= / u £ {x)- ( 9(x, — ) — / 6>(x, v)cfo ) dx 
t^M e\ e J Jd J 

= £ d I u £ (xk + ew)-(d(xk + ew,w) - 9{xk + ew,v)dv 
k% hd £K hd 

= V e d I -({u £ 9)(xk + ew,w) - [ (u £ 9)(xk + sp,w)dp)dw 

JT d e V JT d ' 



We employ the Poincare inequality 



v(z)- / v{z')dz' 



< / \Vv\{z')dz' 



for v(z) = {u £ 9)(xk + sz,v) to obtain 



{u e 9)(xk + ez,v) - I (u £ 9)(xk + £P,v)dp 



< / \V x (u e O)\(xk + £PiV)dp 

JT d 



and 




\V x (u £ 9)\(xk + ep,w) dpdw 



< \\u F 



ijd jjd 

\V x (u £ 9)\(x, v) dxdv 

rllff^fi) (ll^lli 2 (f7xT d ) + II Va:<9||x 2 (OxT d ) y 

From here we obtain (|2.20|) -? . (|2.20p q (along a subsequence), and that 

g e H- 1 (n;L 2 (T d )) 

In addition, we have 

V„5 e = e(V x + -V v )g £ - eV x g £ 

= (V x u £ )(5 p (v - -) — l) - eV x g £ 
e b L 2 (n;M p )+H- 2 (n-M p ) 
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Therefore, 

<V v g £ ,9>= / V x u £ (x)(0{x,-) - / 9(x,v)dv)dx + e < g £ ,V x 9 
Jn \ e J T d / 

| < V v g, > I = lim I < V v g £ , 9 > | 



> 



< lim 



\V x u £ (x)9(x, -)\dx + 



I V x u £ \ \6\dv dx 



< C|l#llL 2 (nxT d ) 

which implies 

V„g G L 2 (Q x T d ) 

and gives using the Poincare inequality the desired (|2.2U|) a. 

To see the first property, consider a test function 9 = 92 <g) ip which is a 
tensor product of ip G C£°(ft) and tp G C°°(T d ). Then 

1 x 

< - -) - 1),<P® V> > 



k£Z d 
k& d 



1 



d [ ¥{Xk + zw) - ip(xk) ~ ^f(Xk) ■ ew 



(ijj(w) — / tp)dw 



ip)dw 



+ ]r e d v^( X fc)- f wfaw)- 

kez d Jd 
0(e)+ / Vip(x)dx ■ / w(ip(w)- / V) 

./o .ITd V ./Td ' 







as ip is of compact support. Since 

l( 5p (v-Z)-i) G b H- 1 (n;M p ) = (H^n- 1 C p ))* 

and finite sums of tensor products YljPj ® V'j are dense in Hq(Q;C p ) we 
obtain (^2D|) i. □ 



3. HOMOGENIZATION OF HYPERBOLIC EQUATIONS 

In this section we consider certain homogenization problems for transport 
equations. First we develop an example where the effective equation can be 
calculated with the help of the double scale kinetic decomposition. Then 
we provide a counter-example where the double scale limit is not the right 
object to treat the effective equation. 
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3.1. Effective equation. Consider the transport equation 

du £ ( x, 

— — + a(-) ■ V x u £ = 

(3.1) dt £ 

u e (0,x)=U°(x,^) 

We assume that a(v) is a C 1 vector field, periodic with period 1, and satis- 
fying diva = 0, and that the initial data U° G L 2 (M. d x T d ) is 1-periodic in 
v and satisfy the uniform bounds 

(hd) f \U°(x,-)\ 2 dx<C. 

Under this hypothesis standard energy estimates for (|3.ip imply the uniform 
bound on solutions 

(3.2) f \u £ (t,x)\ 2 dx< [ \U°(x, -)\ 2 dx < C Vt>0. 
jR d J^d e 

We introduce 

(3.3) fe{t,x,v) = u £ {t,x)5 p (v - j) t x em. d , v eT d , 

and use Lemma 12. II to check that f e satisfies 

df 1 

(3.4) -g + a(v) ■ V x f £ + -a(v) ■ V v f £ = 

(3.5) f £ (t = 0,x,v) = U u (x,v) 5 p (v-j) 

with periodic boundary conditions (v G T d ). The uniform bound (|3.2p 
implies 

(3.6) f £ e b L°°{[0, oo), L 2 (R d ,M p )), 

and thus by Proposition 12.31 along a subsequence if necessary, 

(3.7) f £ ^f weak-* in L°°([0, oo), L 2 (R d ,M p )) 
with / enjoying the improved regularity 

(3.8) / G L°°([0, oo), L 2 (R d x T d )) . 

Our objective is to calculate the effective limit of (|3.ip by computing the 
hydrodynamic limit problem for the kinetic equation (|3.4p - (|3.5p . Note that if 
/ satisfies a well-posed problem then this provides a complete determination 
of the weak limit of u £ since 



fedv ^ fdv = u 
jd Jjd 



We introduce 

(3.9) K = {o G L 2 (T d ) a(v) -V v g = inV'Y 
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and remark that K is the space of solutions of the cell-problem obtained by 
the method of multiscale asymptotic expansion [5] for the homogenization 
problem (|3.1h (see [15], [9]). Let P denote the L 2 -projection operator on the 
kernel K. We prove 



Theorem 3.1. Under hypothesis flh^D the effective limit of problem (I3.ip is 
obtained as u = J Jd fdv where f G L°°([0, oo) , L 2 (R d x T d )), f(t,x, •) G K 
for a.e. (t,x), and f is the unique solution of the kinetic problem 

% + (Pa) ■ V x f = 

(3.10) dt y ' XJ 

f(t = 0,x,v) =PU°(x,v), 
where P is the projection operator on the kernel K . 

Proof. Let f £ and / be as in (|3.6|) , (|3.Tj) and ()3.8[) . The proof is split in three 
steps: 

Step 1 : The limit f belongs to K. The kernel K is defined in (|3.9[) . We 
may consider elements of K as functions of t, x and v instead of only v , as 
t and x play the role of parameter in the definition of K. Thus we have 

(3.11) K x = G L 2 (R d x T d ) a(v) -Vv9 = inV'V 
and 

(3.12) K t . x = |g £ L°°([0, oo), L 2 (M d x T d )) a(v) -V v g = in V'X. 
We may also define all the 

K? = [g G LP(T d )) \a(v)-V v g = 0inV'}, 

and their extensions K x and Kf x . 

The convergence (f37T|) states that for in L^fO, oo), L 2 (M d , C p )) we 
have 

/ / <K*> x, v) df e — > I i (j)(t, x, v) df. 

J0 JM. d xT d JO jR d xT d 

Take (p G C£°([0, oo) x M d x T d ) and compute 



oo 



a(v) ■ V v <f>df £ = -e / (d t 4> + o(u) • V x <f>) df e . 

JK. d xT d 70 Jm d xJ d 

Passing to the limit, we conclude that 

a{v) ■ V v <t>df = 0. 

JR d xT d 

On the other hand / G L°°([0, oo),L 2 (M d x T d )), so / G 
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Step 2 : The limit equation. Consider a function cf> G K x . We wish to 
mollify <f> and use it as a test function in the weak form of (|3.4p . Since 
K x depends only parametrically in x, we may select (f> to be compactly 
supported in x. 

Take H(x) G C c °°(R d ) with f Rd H(x)dx = 1, and fl-(u) G C* c °°((0, 
with J" Td -ff (f ) cfo = 1. For any eft define 

n d H(n(x-y)) <j)(y,v)dy, 
H m (v - rf) 4> n {x,rf) drj, 



with H m {v) = m z2kez d H(m(v + k)), periodic and well defined for all m 
as H(mv) is compactly supported in (0, l/m) d . 
Then for any cj) G K x we have 

/ <j>(y, r])a(rj) ■ ^H m (v - r))H n (x - y) dydi] = . 

Thus, for a Lipshitz continuous, 

a(v) ■ V v 4> n>m df £ 

i x jd 

(a(v) - a(rf)) ■ VH m (v - rj) cp n (x, rj) drj df £ 



l x7 d JT d 



[ [ f 1 (C • Va(« - (1 - t)0) 

JM. d xT d JT d JO 



■ m d+1 VH (m(C + k)) <j) n (x, v - () dt d( df e . 

k&Z d 

Notice that m d+1 (®VH(m() converges in the sense of distributions toward 
C (Id) 5 with C a numerical constant. Moreover thanks to (|3.6p . and to the 
fact that 4> n G L 2 (T d , C c (R d )) and Va G C{T d ), we may pass to the limit 
in m in the previous equality and find 



m— >oo 



Jim / a(v) ■ V v (f) nim df e = C I diva(v) 4> n (x, v) df e = 0, 



as a is divergence free. Multiplying (|3.4p by <j) nm and taking first the limit 
?n — > oo and then the limit e — > 0, we find that for any G compactly 
supported in x we have 



/ <j>nf dxdv - I a(v) ■ V x (j) n f dxdv = 0. 

jR d xT d JR d xT d 

This relation can now easily be extended by approximation to any (ft G K x . 
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Let us denote a the orthogonal projection on K of a. The new function 
a belongs to L°°(T d ) as the projection operator P is continuous on every 
L p (Y d ) for all 1 < p < oo, but does not necessarily have any further reg- 
ularity, Lipschitz for instance (see the appendix where we recall the basic 
properties of P). Now as / G K tjX and V x <f> n G K x , then V x 4>n f £ K} x and 
consequently 

If 

dt jR d xJ d 

On the other hand, the projection operator P may be trivially extended on 
Kt t x from K as t and x are only parameters and of course it commutes with 
derivatives in i or Now, for any G L 2 (R d x T d ) 

/ O - P<f>) f dx dv = 0, / a-V((t) n - P(f) n ) fdxdv = 0. 

7R d xT d ./R d xT d 

Finally for any G L 2 (R d x T d ), we have that 



4> n f dxdv — / a(w) • V x 4>n f dxdv = 0. 

jR d xT d 



dt / 0n/ dxdv - / a(u) • V x </>„ f dxdv = 0. 

This implies that / G is a solution in the sense of distribution to 
(3.13) 8 t f + d(v) ■ V x f = 0. 

Step 3 : Conclusion. Let us begin with the identification of the initial 
value f(t = 0) which has a sense since d t f G L°°([0, oo), H' 1 ^, L 2 (T d ))) 
because of (f3TT3l) and as / G L°°([0, oo), L 2 (M d x T d ))). For every (j) £ K x , 
as 

(3-14) ^ / 4> n df E = f a(v)-V x 4> n df e , 

at JR d xT d JR d xT d 

then f R d x jd 4>ndf £ {t, ., .) has a limit as t — > and this limit is, thanks to 
([33]) 

/ (j) n (x,x/e)U (x,x/e)dx. 

JR d 

Moreover because of (]3.14p . we may pass to the limit in e and deduce that 
(j) n f(t,x,v)dxdv — > / <p n f(0,x,v) = / (p n U°(x,v). 

l xJ d JR d xT d JR d xT d 

On the other hand we of course have for any (j> G L 2 as / G Kf X 



= / ((p n - P(j) n )f(t,X,v)dxdv -> / {<p n ~ P(j)n)f(0,X,v) 

JR d xT d JR d xJ d 

Combining the last two equalities we get that 
(3.15) f{t = 0,x,v) = P U°(x, v). 
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Finally, we notice that Eq. (|3.13p combined with (|3.15p has a unique 
solution in the space of distribution, through standard arguments of kinetic 
theory and as, even though a is only bounded, it does not depend on x. 
Therefore any extracted subsequence of f £ has only one possible limit and the 
whole sequence f £ converges toward the solution of f)3. 13|) with f)3. 15j) . □ 



Examples. We calculate the equation for the double scale limit / and 
the associated effective equation for certain examples, always within the 
framework of (|3.ip , 

1. First consider the case that a(y) is ergodic. Then 

K = {g G L 2 (T d ) : g = const.} 



Pk9 = / gdv =: g 

Jl d 

The equation for / becomes 

d t f + a-V x f = 

and of course u = J Td fdv. 

2. Consider next the homogenization problem 

d t u £ + b(-^)d Xl u £ = 

u £ (0,xi,x 2 ) = U°(x 1 ,x 2 , ■ 

where u e = u £ (t,x), x = (xi,x 2 ) G M 2 , and the vector field a(xi,x 2 ) = 
(b(x 2 ), 0) corresponds to a shear flow with b(x 2 ) ^ for a.e. x 2 . We compute 

K = {g e L 2 (T 2 ) : b(v 2 )d Vl g = 0} = {g = ^(v 2 ) I G L^T 1 )} 



Pk9 = \ g(vi,v 2 )dv 2 
Jo 



Since / G K we conclude that / = f(t,xi,x 2 ,v 2 ) and satisfies the problem 
d t f + b(v 2 )d x J = 



f(0,xuX2,v 2 ) = P K U°= I V 

Jo 



(xi,x 2 ,v 1 ,v 2 )dvi 



The weak limit u = J J2 f satisfies the integrated equation. 

3. It is possible to give a more general framework for the situation of the 
previous example. Suppose that the divergence free vector field a is such 
that the following description of K is true: There exist functions £l, . . . (,n, 
N < d from to R. These functions are local coordinates in the sense that 
they may be completed by £/v+i> • • • £<2 and that the change of coordinates v 
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to {£i(v), . . . , £d(v)) is a C 1 diffeomorphism from T d to some domain O C R d . 
And finally 

K = MSiiv), . . . , £ N (v)) | G ^ 2 (0)}. 

For instance in dimension d = 2, as div a = 0, there is always £ : M d — > R 
such that a = V" 1 "^. Now if in addition £ is a periodic regular function with 
V£(v) / for all v, which is a non trivial assumption, then K is exactly the 
set of functions 

In that case, we may define g(t, x, £) = f(t, x, F(£l, . . . , £<z)) with 1/ the 
inverse change of variables. Then g does not depend on £jv+i, ■ ■ ■ , £<z and it 
simply satisfies 

% + &(&,..., frr)-V x = O, 

with 6(£) = o(y(0) = / a(F(£))^ A r + i . . . d£ d . 

4. Notice now that the kernel K, endowed with the usual L 2 scalar 
product, is a Hilbert space and so that the kernel K admits an orhonormal 
basis {i>k(v)}, possibly countable. Since / G K it will be given in a Fourier 
expansion 

oo 

/ = X m fc(*> zOV'fcfa ) where m k =< /, ip k >. 
k=i 

Moreover, we see that 

< P K {a)f,^ k >=< P K {af),ip k >=< af,ip k > 

and one computes that the set of moments m k satisfies the initial value 
problem 

d cv q 

d t m k + X ( X < a i^ n > ^ k > = 

j=l n=l 3 

m k (0, x) =< P K U°(x, ■), i/> k >= / U°(x, v)ip k (v)dv 

As the wave speed a is real, < ajip n ,^k >=< ip n , a-j^k >, and the system of 
moments is an infinite symmetric hyperbolic system. 

3.2. The multiscale case: A counterexample. A natural extension of 
the previous analysis is to deal with transport coefficients depending on more 
than one scale. Consider for example the equation 



(3.16) 



d t u £ + a £ ■ V x u £ = 0, 
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with a £ = a(x,x/e) and a(x,v) a Lipschitz function, or even with a £ = 
a(x,x/e,x/e 2 ) (or with as many scales as one cares to introduce). Assume 
again that div v a(x,v) = and div w a(x,v,w) = 0. 

Is it possible to derive an equation for the double scale limit (or for the 
triple scale limit when a(x,x/e,x/e 2 )) in the case of (|3. 16|) ? 

In fact, it is relatively easy to show that the previous approach does not 
work! Everything goes as before in the beginning; upon defining 

f e (t,x,v,w) = u £ (t,x) 5 p (v - x/e) 5 p (w - v/e), 

as in paragraph 12.1.31 one simply obtains the generalized kinetic equation 

dtfe + • (a(x, v, w)f e ) + -a ■ V v f £ + \a ■ V w f £ = 0. 

However it is not always possible to derive a well posed problem for the 
hydrodynamic limit, even in the simple setting where a depends only on x, v 
and the equation for f e (t,x,v) is 

dtfe + V x • (a(x, v)f £ ) + -a(x, v) ■ V v f £ = 0. 

Indeed the only information that we have is that any limit / belongs to the 
kernel which now depends on v and x 

K = {f e L 2 (R d x T d ) | a(x, v) • V v f = 0}. 

On the other hand, when projecting the equation on K, it is not possible 
to handle the term with the x derivative as projection on K and differenti- 
ation in x no longer commute. This is associated to the possibility that the 
dimensionality of K may vary with x. 

In addition, at the level of the double scale limit, this is not a mere 
technical problem, rather the double scale limit is in general not unique and 
depends on the choice of the extracted subsequence in e. 

This can be simply seen for the problem 

d t u £ + a(x) ■ V x u £ = 0, t £ E + , x G R 2 , 

(3.17) 

u £ {t = 0,x) = U°{x,x/e). 
The oscillations are due only to the initial data as the transport coefficient no 
longer depends on e, and u the weak limit of u £ satisfies the same equation. 
Take now 

ai{x) = 1, a 2 {x) = x x \< X1 <\ + 1 1W>1> 
so that a is Lipschitz and divergence free, and select the initial data 

U°(x,v) =K(x)L(v 2 ), 
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with K and L two C°° functions, L periodic on R of period 1 and with zero 
average, and K compactly supported with support in x\ in [— 1, —1/2]. 

As the average of U°(x,v) in v vanishes for all x, the weak limit u of u £ 
is uniformly 0. 

For 1 < t < 3/2, the support in x\ of the solution u £ is entirely in the 
interval [0, 1]. Therefore any double scale limit / should satisfy 

d vi f + x 1 d V2 f = 0. 

It is easy to check that the only L 2 solutions to this last equations are the 
functions which depend only on x and not on v. Therefore for 1 < t < 3/2 
the double scale limit is equal to u, i.e. uniformly vanishes. 

Let us finally compute the double scale limit for t > 2 and check that it 
does not vanish. For that introduce the characteristics X(t,x) 

d t X(t, x) = a{X(t, x)), X(0, x) = x. 

We need the characteristics only for those x which belong to the support of 
U°, that is for -1 < x± < -1/2. 
As a\ = 1, we simply have 

Xi(t, x) = x\ + t. 

As to X2, as long as X\ < or t < to = —x\ (remember x\ G [—1, —1/2]) 
it is equal to X2- For to < t < t\ = 1 — x\ (corresponding to X\ in [0, 1]), 
we have 

d t X 2 =X 1 =x 1 + t. 

As a consequence 

t 2 t 2 (l-Xi) 2 x 2 
X 2 (h) = x 2 + Xl (h - t ) + -j - -j = x 2 + Xi + ^ -^=x 2 + 1/2. 

After t > ti, dtX2 = 1 and so 

X 2 (t) = x 2 + 1/2 + t - h = x 2 + 1/2 + t - 1 + xi = x 2 + X 1 - 1/2. 

With this, the solution u £ is given for t > 2 by 
u £ (t, x) = u £ (0, xi — t, X2 — Xl + 1/2) 

= K(xi -t,x 2 -xi + 1/2) L(x 2 /e - x 1 /e + l/2e). 

For every a £ [0, 1], choose a subsequence e n such that l/2e n converges to 
a modulo 1. Then the double scale limit associated to this subsequence is 
the function 

f(t, x, v) = K(xi - t, x 2 — xi + 1/2) L(v 2 - v\ + a). 
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Instead of one unique limit, we obtain a whole family which clearly indicates 
the ill-posedness of the problem at the level of the double scale limit. 

4. Transport of oscillating fine-scale 

An interesting question that can be studied using the techniques devel- 
opped in section [3] is the problem of transport of an oscillatory fine-scale 
structure under a divergence-free vector field. Consider the homogenization 
problem 

d t u £ + V x ■ A(t, x, ^^-)u £ = 

(4.1) e 

u £ (0,x) = U°(x, ^) 

where ip : R d x R — > R rf is a C 2 map describing the fine scale of oscillations 
that satisfies for some c > 

(p(-,t) is surjective and invertible for t fixed 

(h so ) det VO/T 1 (-,£)) >c>0 

ip(-, 0) = id 

and a e = A(t,x, is a divergence free field. 

The latter is guaranteed provided A(t, x, v) is a C 1 vector field 1-periodic 
in v such that 

V x -A(t,x,v) =0 

tr(V^V^) = ^— — = 
i,j=l 

The initial data [7° G L 2 (R d x T d ) are 1-periodic in v and satisfy the uniform 
bounds 

(hd) / |t/ (x,-)| 2 dx<C. 

Then standard energy estimates imply that solutions of (|4.ip satisfy the 
uniform bound 

(4.2) / \u £ (t,x)\ 2 dx < [ \U°(x,-)\ 2 dx < C. 

JR d jR d £ 

Our objective is to calculate an effective equation for the weak limits of {u £ }. 
The counterexample of section 13.21 indicates that we can not expect to do 
that in full generality. A more precise statement of what will be achieved 
is that we will identify conditions on the vector field A and the structure 
function ip under which an effective equation is calculated. 
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4.1. Reformulation via a kinetic problem. We introduce the "kinetic 
function" 

(4.3) fs = U£ (t,x)5 p {v-?^-), 

which is well defined (see section 12. 1 ,2[) as a measure. Due to the identities 
in lemma 12.21 of section 12.1.21 it is possible to transform the homogenization 
problem (|4.1|) into a hyperbolic limit for a kinetic initial value problem: 

Lemma 4.1. If u £ a weak solution of (|4.ip then f £ in (|4.3p verifies in T>' 
the kinetic problem 

dtf £ + V x -(Af £ )+V v (-Bf £ )=0, 

(4.4) £ 

f £ (0,x,v)=U°(x,v)5 p (v-^), 

where the vector field B(t,x,v), defined by 

(4.5) Bi = [d t + A- V x )<f>i, i = l,...,d, 
is 1-periodic in v and divergence- free, V v ■ B = 0. 



Proof. Let u £ be a weak solution of (|4.ip . By (h tv f), B defined in (|4.5p 
is 1-periodic in v and divergence free. We consider a test function 8 G 
C£([0,oo) x R d ;C 1 (T d )) and compute 

f°° f f [d t 9 + V x (A-6) + -V v (B-9)]df £ 

Jo JR d Jf d £ 

= f f u £ (x,t)(d t 6 + A-X7 X 8 + -B ■X7 v 9)(t,x,-)dxdt 
Jo JR d v e J e 

= ^ [ u £ \dt{9(t,x,£)) +A(t,x,£)-V x (0(t,x^)) 
Jo JR d Lee e 



U (x,-)e(0,x,-)dx 

d £ £ 



dxdt 



8(0,x,v)df £ (0,x,v); 

R d Jj d 

that is, f £ is a weak solution of (|4.4p , □ 

An alternative, albeit formal, derivation of (|4.4p may be obtained by 
studying characteristics. The characteristic curve of (|4.ip emanating from 
the point y is defined by 

1 a;(0,y) = y 



26 



P.-E. JABIN AND A.E. TZAVARAS 



and is denoted by x = X(t; y). Along such curves we have 

=i(^ + A(*, t> f).V.* 

The two equations together can be embedded into the system of ordinary 
differential equations 



(4.6) 



f = i (ip t + A(x,t,v)-V x ip 



in the following sense: If (Y(t;y,u),U(t;y,u)) is the solution of (|4.6|) ema- 
nating from the point (y, u) then 

X(t;y)=Y(t;y,^(y,0)) 
<f{X(t;y),t) 

= U (t; y, <p(y, 0)) 

Note that (|4.4|) i is precisely the Liouville equation associated to the charac- 
teristic system (I4.6p . 



4.2. Conditions leading to an effective equation. Our next goal is 
to derive an effective equation for the hydrodynamic limit of (|4.4I) . We 
first show that under hypothesis flh so D the definition (|4.3p still induces good 
properties for the weak limit points of {/ e }. 



Lemma 4.2. Under hypotheses flh so D and (|4. 2 j) . 

(4.7) Sp (v-^^)^1 inV' 

(4.8) f £ G b L°°((0,^);L 2 (R d ;M p )) 
and, along a subsequence (if necessary), 

(4.9) f e -^f weak-* in L°°((0, oo); L 2 (R d ; M p )) 

(4.10) with f G L°°((0, oo); L 2 (R d x T d )) 

Proo/. For G C c °°((0, oo) x R d ; C°°(T d )) we have 

<S p (v-^-),e>= r [ 9(x,t,^-)dxdt 
£ Jo Jm d £ 

8{cp- 1 {y,t),t,V-)\detX7 v (<p- 1 )\dydt 

£ 



9(t, x, v)dvdxdt 

Note next that for (x, t) fixed 

\\f £ (x,t,-)\\ Mp = \u £ (x,t)\ 
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and thus (gSJ) and gJl follow from Finally, 

I < / £ ,0 > | < ||w e ||i«(i,2) ^ (J^\9(x,t,^^-)\ 2 dx^ dt 

I (I I \0(x,t,v)\ 2 dxdv) dt 

JO \JR d JT d J 

and KW\i follows. □ 



Remark 4.3. A hypothesis of the type of (\h. so \j is essential for the validity of 
(|4.7p and accordingly for (|4.10p . For instance, in the extreme case that ip is 
a constant map, <p(x, t) = c, it is possible by choosing appropriate sequences 
e n — > to achieve any weak limit 

5 p (y ) — 5 p (v — v ) with any < v < 1. 

The regularity of / is then no better than the regularity of {f e } and (|4.10l) 
is of course violated. 

We conclude by providing a formal derivation of an effective equation. 
Consider the e — > limit of (|4.4p -( l4~5j) and recall that, by (14. 2p and lemma 
B~2l we have f e / as in (fOj) and ([4TTD]) . Define the set 

Ki,* = {S6 L°°((0,oo);L 2 (M d x T d )) | fl(t,a;,v) • V v g = 0} . 

The set K = N(B) is the null space of the operator B := B ■ V„ and in 
general it will depend on (t, x). We will derive the effective equation under 
the hypothesis 

(H) N(B) is independent of (x, t) 

Then we have the decomposition 

L 2 = N(B) © R(B T ) =K®K ± 

and the spaces remain the same for any point (x,t). Let P denote the 
L 2 -projection on the set K. Any 6 6 L 2 (T d ) can be decomposed as 

9 = P9 + (I - P)9 =: ij) + <f> 

Moreover the differentiation operators dt and V x commute with the projec- 
tor P. 

For V G AT (5), using V x • A = V„ • -B = 0, we derive from that 
(4.11) dt<f,ifi> +V X < A/, ^ >= 
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where the brackets denote the usual inner product in L 2 (T d ). One easily 
sees that f(t,x, •) 6 K for a.e. (t,x). Given 9 E L 2 (T d ) let ip = P9. Then 

</,0>=</,^> 

and 

< Aif,ip >=< p(Af)^>=< p(Af),e> 

Since f £ K we have P(A{f) = P(Ai)f and we conclude that (|4.11j) can be 
expressed in the form 

d t <f,9> +V X - < P(A)f, 9>=0, 9 £ L 2 (T d ) . 

The effective equation thus takes the form 

(4.12) d t f + V x -(PA)f = 0. 

The above derivation of (|4.12p is formal and is based on hypothesis (|H|) . 



which is quite restrictive especially when viewed together with (h tv f) that 
has to be satisfied simultaneously. We view this equation as a theoretical 
framework of when an effective equation can be computed. To derive it 
rigorously one needs an analysis as in Theorem 13 . 1 1 and we will not pursue the 
details here. The counterexample in section [3721 indicates that the hypothesis 
(}H|) is essential. 

We list two examples that can be viewed under the above framework. 
First, the homogenization problem (|3.ip is a special case of (|4.ip with the 
obvious identifications. A second example is given by the problem 

d t u £ + a(x) ■ V x u £ = 

^ ^ u e (x,0) = U(x, ^) 

where a is a divergence free field, V x ■ a = 0. Define ip(t, x) to be the 
backward characteristic emanating from the point x. Then ip = (ipi, ip n ) 
satisfies 

d t <fi + a(x) ■ V x pi = 
ipi(0,x) = Xi 

The problem (|4.13p fits under the framework of (|4.ip under the selections 

A(t, x, v) = a(x) . P>i(t, x, v) = (dt + a(x) ■ V x )(fi = . 
The kinetic equation for f £ = u £ 5 p (y — becomes 

dtfe + a{x) ■ V x f £ = 
f £ (0,x,v) = U(x,v)5 p (v - j) 
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while the limiting / satisfies the same transport equation with initial con- 
dition f(0,x,v) = U(x,v). Hence, it is computed explicitly by 

f(t,x,v) = U(cp(t,x),v) . 

5. Enhanced diffusion 

In this section we study the enhanced diffusion problem 

i x 

d t u £ + -a(x, -) • V x u £ = aA x u £ 

(5.1) £ £ 

u £ (0,x) = U°(x,^) 

where a(x, v) is a Lipshitz vector field periodic (with period 1) in v that 
satisfies 

(h. v f) V x ■ a = V v ■ a = , / a(x, v)dv = , 

a > is constant and U° G L 2 (R d x T d ). We use this as an example to 
develop the methodology of section [2j For previous work and a commen- 
tary on the significance of this problem we refer to Avellaneda-Majda [1], 
Fannjiang-Papanicolaou [10] and references therein. It is assumed that the 
initial data oscillates at the scale e and satisfy the uniform bound 

(hd) [ \U°(x,-)\ 2 dx<C. 

Standard energy estimates then imply the uniform bounds 

(5.2) f \u £ {t,x)\ 2 dx + a [ [ \Vu £ \ 2 dxdt < [ \U°{x, -)\ 2 dx < C 

jR d JO jR d jR d £ 



for solutions of (15. lj) . 

We introduce the kinetic decomposition 

(5.3) f £ (t,x,v) = u £ (t,x)5 p (v - t £ M+, x £ R d , v £ T d , 

and use Lemma 12.11 to obtain that f £ satisfies the transport-diffusion equa- 
tion 

+ -a(x, v) • V x f £ + ^ { a (x, v) ■ V v f £ - aA v f £ ) 

(5-4) = aAx f £ + — V x ■ V v f £ , in V 

e 

f £ (t = 0,x,v) = U°(x,v) 8 p (v-^) 

with periodic boundary conditions on the torus in the v variable. Our ob- 
jective is to analyze the e — > limit of this problem and through this process 
to calculate the effective equation satisfied by the weak limit of u £ . We note 
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that this is a hydro dynamic limit problem in the diffusive scaling for the 
kinetic equation (|5.4p . 
We prove. 



Theorem 5.1. Under hypothesis (h v f ) and (h^) we have the following as- 
ymptotic behavior for f e as e — ► 0: 

(5.5) f e (t,x,v) ^ u(t,x) weak-* in L°°([0, oo), L 2 (R d ,M p ) 

(5.6) f E (t,x, v) = u e (t, x) + eg(t, x, v) + o(e) , m V' 

where u of class (|5.2p and o G £ 2 ((0, oo) x Q; ff 1 (T d )) satisfy respectively 

(5.7) <9tii — aAjii + V x ■ / a(x,v)g(t, x,v) = 

(5.8) ctA v g — V v • ag = a ■ V x u 
The weak limit u satisfies the effective diffusion equation 

d t u = a d x . ^(<% + ^ V„Xi • ^vXjdvjd x 



h3 



u(0,x) = / J7°(x,u)cfo 
where Xk> k = 1, ...,d, is the solution of the cell problem 
(5.9) aA v Xk -V v -axk = a-e k 

Proof. Let f E be defined as in (|5.3p . Then / e satisfies the problem (|5.4|) and 
"Ue = Jjd fe- The proof is split in three steps: 

Step 1 : Characterization of the weak limit. From (|5.2p and Lemma 12.11 
we obtain uniform bounds for f e : 

f E e b L°°([0, oo); L 2 (R d ,M p ), 

( 5 - 10 ) 1 

(V, + -V v )f £ e b L 2 ((0,oo) x R d ; M p ) 

Using (a slight variant of) Proposition 12. 31 we see that, along a subsequence 
if necessary, f e satisfies 

f e -^f weak-* in L°°([0, oo); L 2 (R d ,M p ) 

[>A ' * fe L°°([0, oo) ; L 2 (R d x T d )) , 

that is 

J 6(t, x,v)df £ (t, x,v) — » y 9(t,x,v)f(t,x,v)dtdxdv 
for (9 € L 1 ((0,oo);L 2 (R d ,C , p )). 
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Passing to the limit e — > in (15. 4p and using (15.101) we see that / satisfies 

aA v f -V v -af = 0, inV, 

and for any test function 8 

<V v f,9>= lim<V v f £ ,9>=0. 

£— >0 

Hence, V„/ = in T>' and 

(5.12) f(t,x,v) = [ fdv=:u(t,x) 



Step 2 : Asymptotics of f £ . Define next 

(5.13) g £ (t,x,v) = -(f E (t,x,v)-u £ (t,x)), t G K+ , x G R d , v G T d , 

where u £ = J" Td / e . We proceed along the lines of Proposition 12.51 replacing 
the bounds of that proposition by the bound (|5.2|) and accounting for the 
extra dependence in time. After minor modifications in the proof we obtain 
for any T > 

g £ G b L 2 ((0,T); H-H^Mp)) 

( 5 . 14 ) 9s^9 weak-*inL 2 ((0,T); H^(n,M p )) 

5 GL 2 ((0 iro )xfi;^(T d )), / 5 = 0. 
Accordingly, f £ enjoys the asymptotic expansion 

f £ = u £ + eg + o(e) in X>' . 



On the other hand, on account of (15. lj) . (|5.4p and (ht,j), it follows that 
u £ and g £ satisfy 

/• f £ — u £ 

(5.15) d t u £ - aA x u £ + V x - j a(x,v)— - = 

and 

(5.16) aA v g £ - V v • ag £ = e(d t f £ - aA x f £ ) - 2aV x • V v f £ + V x • a/ e 

Using (|5.1ip . (|5.12p and (|5.14p . we pass to the limit e — > and deduce that 
it, <? satisfy (|5.7p and ()5.8[) respectively. 

5fep 3 : Characterization of the limit problem. Due to its regularity the 
solution g of (15. 8p is unique and can be expressed in the form 

g = V x u(t,x) • x(x,v) 

where \ = (Xli •••>Xd) is the solution of the cell problem 

(5.17) aA v Xk -V v -axk = a-e k = a k . 
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A direct computation shows that solutions of (15. 9b satisfy the property 

/ ( a iXk + a k Xi)dv = -a [ V v Xi ■ ^vXkdv 

1 Jjd Jjd 

and (|5.7p may be written in the equivalent forms 



hi 
a 



The latter is a diffusion equation with positive definite diffusion matrix 

Dij = Sij + / V v Xi • ^vXjdv 
J2D ij ^ j = \u\ 2 + [ \V v { X -v)\ 2 dv, ueM. d , 

determined through the solution of (|5.9|) . □ 

6. Appendix I 

We prove in the appendix a lemma that is used in the justification of 
multiscale decompositions. 

Lemma 6.1. Let Q be an open subset ofR d , 9 £ C C (Q), ip G C(T d ), if) G 
C(T d ), and suppose that 5 = 5(e) —* as e — > 0. Then, as e — > 0, 

(6.1) / 9(x)ip(-)dx^ [ <p(z)dz [ 6(x)dx 
Jn £ Jf d Jn 

(6.2) / e(x)(p(-)ipAdx -> / V(y)^y / f 9(x)dx 
Jn e eo Jr d Jj d Jn 

Proof. Fix 9 G C c (fl) and let e < dist (supp6>, <9f2). We consider a cover 
of the support of the function 9 by cubes Ck centered at points \k G 
of latteral size e. The number of the cubes covering supp 9 satisfies Ns d = 
O(|supp0|), and since e\fd < dist (supp#, d£l) the covering can be arranged 
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so that supp6> C uf =1 C fc C U while \xk G %> d - We have 

r N r 

/ 6(x)ip(-)dx = J2 G(x)ip(-)dx 

TV 1 

= YV / 0(xfc + £*) y(-Xfc + *) d2 



k=l 
N 



k=i 

and thus 

I = 6{x)ip{—)dx — / <p(y)dy / 6{x)dx 
Jn £ Jf d Jn 

= J2 £d e( Xk + ey)if(y)dy - ip(y)dy YV / ff(x fc + e«)(fa 
fc=i ^ Td ^ T<i fc=i 

TV 

= / <p(v)J2e? {0{xk + ey)-6{ Xk + ez))dzdy 

JT d k=1 JJ d 

Using the uniform continuity of 6 and that N = 0(-^) we deduce I — > as 
e -> and (loTTD . 
Next observe that 

0(x)ip(— )ifj(-^-)dx — / ifj(y)dy / (p(z)dz / 9{x)dx 



X X i 

(x)<p(-)ip(— )dx - \ tp(y)dy \ 9(x)ip(-)dx 
e ed J T d 



+ 

/1 + /2 



4>(y)dy ( / 6(x)(p(—)dx — / tp(z)dz / 6{x)dx 
i d \Jn e JT d Jn 



and that 72 — ► as e — > 0. 

Consider now a covering of supp# by cubes C& centered at points Xfc £ 
e<5Z d of latteral size e5. As in the preceding argument we can arrange the 
cubes so that supp# C U k=l Ck C fi and their number TV satisfies N(e5) d = 
0(|supp0|). We have 

/ 0(xM-)^)dx = Y, e(xM-)i;(-)dx 
Jn e e5 ^Jc h £ e6 

N 1 

= V(e<5) d / 9{xk + e5z)ip(-Xk + 6z)ijj(z)dz 
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and 



h 



N N 



,X 



N 

Jf d Jf d v e 



Again, since iV = 0(-r~xy[) an d lim £ ^o <5 = 0, we deduce 1% — ► as e — > 



Remark 6.2. Relation (|6.ip is classical (see [5]) and is only proved here as 
a precursor to the proof of (|6.2p . Equation (|6.2p indicates that oscillations 
of different scales do not correlate and suggests that the definition (|2.8|) is 
meaningful. 

Both equations can be extended for test functions 9 £ C(0) provided 
that Q C R d is a bounded open set and its boundary dfl is of finite d — 1 
Hausdorff dimension. 

As an application we prove (|2,14p . 

Lemma 6.3. We have 



Equation (|6,ip justifies that for 9 = x( x ) ® f{ v ) a tensor product with 
X S C C (J7) and 99 £ C p and by a density argument also for x £ 

To complete the proof we need to show that finite sums of tensor prod- 
ucts YljXj ® y>j are dense in L l (Q,,C p ). Fix G L l (Q,,C p ) and consider a 
decomposition of the torus T d into squares of size 1/n. Take a partition of 
unity </?j S C p , i = 1, ...,n d , with each </?j supported in a square of size 2/n 
and £^ = 1. Let V{ be the center of each square and define 



and ([BT2]) . 



□ 



(5 p (w ) ^ 1 weak-* in L°°(^, M p ) 

Proof. We need to show that for 9 6 L 1 (r2, C p ) we have 
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Clearly, 9 n is a sum of tensor products. Now define 

sup \9(x, v) — 9 n (x, v)\ < sup \9(x,v) — 9(x,v + h)\ =: g n (x) 

vdT d t , 6 T<*,|A[<| 

and thus 

¥ - Qn\\v-(c p ) < I \g n {x)\dx 

J S2 

Note that g n {x) — ► for a.e. x £ £1 and that |p n (a;)| < 2sup 1 , eT d \9(x,v)\. 
The latter is an L 1 function by the very definition of 9, and the dominated 
convergence theorem implies f n \g n \dx — > 0. □ 

7. Appendix II: Some basic results of ergodic theory 

The purpose of this appendix is to recall some well known properties of 
the classical ergodic theory for the projection on the kernel 

K = {/ G L 2 (T d )\ a(v) -V v f(v) = 0}, 

where the last equality is of course in the sense of distributions. 

Let us define the characteristics associated with a which are the solutions 
on T d of the following differential equation 

d t T(t,v) = a(T(t,v)), T(0,v)=v. 

Then assuming that 

(7.1) a€W 1 ' 0O (T i ), V v -a = 0, 

the characteristics are well defined and for a fixed t, the transform v — > 
T(t, v) is a mesure preserving homeomorphism of T d . We then have the 
well-known theorem (see Sinai [IB] for more details) 

Theorem 7.1. For every f £ L p (T d ) with 1 < p < oo, there exists a unique 
function in L p (T d ), denoted by Pf, such that 
1 

7 

Moreover Pf satisfies in the sense of distribution 

a{v) ■ V v Pf(v) = 0, 

and if f E L 2 (T d ), then Pf is exactly the orthogonal projection of f on K. 

This immediately implies the 

Corollary 7.2. The orthogonal projection on K may be extended as an 
operator on L p (T d ) for every 1 < p < oo. In addition if f E L 2 n L p (T d ) 
with 1 < p < oo (p = oo allowed), then Pxf also belongs to L 2 n L p (T d ). 



t 

f(T(s,v))ds — >Pf{v), ast^oo, strongly in L p (T d ). 



o 
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Proof of Theorem 17. 1L This proof exactly corresponds to the one in 

in the particular case which we consider. 
Notice that if, for / G L p (T d ), \ JjJ f(T(s,v)) ds converges to Pf then 
trivially 

Pf(T(t,v))=Pf(v) Vt. 
Therefore we automatically have in the sense of distribution that 

a(v) ■ V v Pf = 0. 

Take now / in LP and assume first that there exists g G LP with a ■ V v g = 
and h G W l > p (T d ) such that 

f = g + a- V v h. 
Then notice that in the sense of distribution 

d t (g(T(t, v))) = a(T(t, v)) ■ V v g(T(t, v)) = 0, 

and so 

g(T(t,v))=g(T(0,v))=g(v). 

On the other hand 

a(T(t,v)) • V v h(T(t,v)) = d t (h(T(t,v))) , 

and therefore 

f(T{s,v))ds = g{v) + 



I ./o t 
i r* 



Consequently in this case j J f(T(s,v))ds converges to g which is unique 
as a consequence. This proves the theorem on the set 

L p = {g + a(v) • V v h{v) | h G W 1 ' 00 ^), g G L p {T d ) with a-V v g = 0}. 

Let us first prove that L2 is dense in L 2 (T d ). If Li is not dense, then there 
exists / G L 2 \{0} orthogonal to L 2 . This implies that for all h G W 1,00 (T d ) 

f(v)a(v) ■ V v h(v)dv = 0, 

or in other words / G K. But K C L2 and / should consequently be 
orthogonal to K, which is impossible. Notice that Pf necessarily is the 
orthogonal projection on K as a ■ V/i belongs to K^. 

Now for any / G LP. If 1 < p < 2, take g n + a-V v h n = f n G L2 converging 
toward / in LP (first take f n G L? and then select f n by diagonal extraction) . 
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We have that 

\\j f f(T(s, v)) ds-^J^ f(T(s, v)) ds\\ LP 

< 2||/-/„||lp + (j + ^j HM^- 

So the sequence \ Jq f(T(s,v)) ds is of Cauchy in L p and hence converges 
to a unique limit Pf. 

Finally if / € L°°(T d ), then / G L 2 (T d ) and \ f* f(T(s,v)) ds converges 
to Pf in L 2 . As the first quantity is uniformly bounded in L°°, Pf G L°° and 
the convergence holds in every L p , p < 00. By interpolation, one eventually 
obtains the desired result for / G L p (T d ). □ 
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